Resonant modes in strain-induced graphene superlattices 
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We study tunneling across a strain-induced superlattice in graphene. In studying the effect of 
appUed strain on the low-lying Dirac-like spectrum, both a shift of the Dirac points in reciprocal 
space, and a deformation of the Dirac cones is explicitly considered. The latter corresponds to 
an anisotropic, possibly non-uniform, Fermi velocity. Along with the modes with unit transmission 
usually found across a single barrier, we analytically find additional resonant modes when considering 
a periodic structure of several strain-induced barriers. We also study the band-like spectrum of 
bound states, as a function of conserved energy and transverse momentum. Such a strain-induced 
superlattice may thus effectively work as a mode filter for transport in graphene. 

PACS numbers: 81.05.uo, 72.80.Vp, 85.30.Mn 
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Graphene is a single layer of carbon atoms in the sp^ 
hybridization state, arranged according to a honeycomb 
lattice^'-. Transport properties in graphene are largely 
determined by its reduced dimensionality, which char- 
acterizes its remarkable electronic properties^. These 
include low-energy quasiparticles with a Dirac-like spec- 
trum and a linearly vanishing density of states (DOS) 
at the Fermi level. Evidence of such an unconventional 
behaviour is to be found in several electronic properties, 
such as Klein tunneling^rl, the optical conductivityiSrJ^, 
and the plasmon dispersion relationi^""— . These have 
been predicted to depend quite generally on applied 
strairiii, following the earlier suggestion that suitably 
deformed graphene sheets could be engineered into nan- 
odevices with the desired electron propertiea^^. For in- 
stance, it has been recently demonstrated that the elec- 
trical properties of epitaxial graphene on SiC strongly 
depend on the local strain induced in graphene by the 
substratei^. One thus expects that a suitable pattern 
of periodically repeating stripes, with alternating values 
of strain, i.e. a strain-induced superlattice, may pro- 
duce coherent effects on single particle transport, de- 
pending on the energy and momentum of the incident 
electrons. Here, we therefore study the possible occur- 
rence of resonant states within a strain-induced super- 
lattice in graphene. 

We consider quasiparticle transmission across N iden- 
tical barriers, each of width ^, the inter-barrier separation 
being also i, such that 2Ni — D (Fig. [T]). Let x denote 
the coordinate orthogonal to the barriers, forming an an- 
gle with the graphene zig-zag direction. Thus, = 
{resp., = tt/2) will refer to a superlattice oriented along 
the zig-zag {resp., armchair) direction. Such a superlat- 
tice is usually obtained via a step-wise varying gate po- 
tential U{x) = U±, with U{x) = U- within each barrier 
[2{m-l)£<x < {2m-l)i,m = 1, . . . N], a.ndU{x) = U+ 
between two neighboring barriers [{2m — 1)£ < x < 2mi, 
m = 1,...A'^]. Here, we will additionally consider a 
nonuniform profile of uniaxial strain e = £(2;) applied 
along the direction, with strain modulus alternating 
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FIG. 1: Schematic plot of the superlattice of A'^ identical bar- 
riers, with £ denoting both each barrier's width and the inter- 
barrier separation, while D = 2NI. Subscript — refers to the 
region within a barrier (labelled II), while subscript -I- refers 
to the interbarrier region (labelled I and III). 



between the values e{x) = e± inside and outside a bar- 
rier, as above. Such a dependence approximates a smooth 
periodic strain wave with period 2£, as a train of sharp 
steps. 

Within each barrier, strain is described by the tensor 
£ = ^e[{l — i')-\- {l-\-iy)A{9)], where e = e{x) is the strain 
modulus, J/ = 0.14 is Poisson's ratio for graphene^*^, and 
A{9) = cos(26')cr2 -\- sin{20)(Tx, with ai {i = x,y,z) de- 
noting the Pauli matrices. Within a quite general tight- 
binding approach'^, strain then enters the electronic prop- 
erties through the dependence of the hopping parameters 
on the lattice vectors^!. Expanding such a tight-binding 
Hamiltonian to linear order in the strain modulus, one 
finds that the low-lying spectrum can still be described 
by a Dirac-Uke Hamiltonian, but now (i) applied strain 
shifts the location of the Dirac points in reciprocal space 
with respect to ±K at the vertices of the first Brillouin 
zone, and (ii) it induces a deformation of the Dirac cones, 
which can be accounted in terms of an anisotropic Fermi 
velocity vp. Specifically, one finds for the Hamiltonian 
under applied strain 



H = hvFU''{0)a-qU{0), 



(1) 



where q = (91,92)^ measures the wave vector displace- 
ment from the shifted Dirac points q^a = ±(ko£(1 + 
i/)cos(26i),-Ko£(l + J^)sin(26'))^, a^ = (1 - X^eja^ {i = 



1,2) take into account of the strain-induced deforma- 
tion of the Fermi vefocity, with A^; = 2k, Aj, = —2k,v, 
U{9) = diag(l,e~*^) is the unitary matrix performing 
a rotation mapping the zig zag direction onto the the 
direction x of applied strain, kq = {a/2t)\dt/ da\ « 1.6 
is related to the logarithmic derivative of the nearest- 
neighbor hopping parameter t with respect to the lattice 
parameter a at zero strain, and k = kq — ^ (cf. Ref. |22|). 
Since the strain superlattice is uniform along the coor- 
dinate orthogonal to the direction of applied strain, say y, 
stationary eigenmodes will be characterized by constant 
energy E and transverse wave vector ky. The station- 
ary Dirac equation associated to Eq. ([1]), with appro- 
priate matching conditions for the quasiparticle spinor 
due to the continuity of its associated current density at 
the barriers' edges, can then be equivalently recast using 
the transfer matrix formalis m^^i^^ . Following Ref. llTl . 
for the transfer matrix across the first, say, barrier in 

i)(2£,0) == e*«Di('^++'^-KM(i), where 
— 1), and M'^' is a unimodular matrix. 
Specifically, one obtains 



Fig. [TJ one finds 
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A + IT], 

s\nh.{q-€) sinh((7_|_£) 



(2a) 



(k_K-|_ — U-U+) 



q- 94 

-cosh(g_^)cosh(g+^), (2b) 

U+U_ - K+K_ . 

smh((7_£) cosh(g+£) 

q+q- 

sinh(q+^)cosh(g_^)j, (2c) 



where A is always real, whereas ry can be real or purely 
imaginary, depending ky and E. More compactly, one 
also finds 



E« 



exp(g+£) 



q+q- 



■ sinh(q_£) + cosh(g_^) 



(3) 
In Eqs. ([2]) and ([3]), we have employed the definitions 

K± = (1 - \£±){ky - q^'ye±)l{l - \x£±), u± ^ {E - 

U±)/[hvF{^ — Ax^i)], and q± ~ \/^±— ""i- Making 
use of the Chebyshev identity for the A^th power of a 
unimodular matrix^^, for the evolution matrix across N 
identical barriers, one finds^^ 

r~i-i>-,w smh(Nz) ~ (-i\ smh((N — l)z) , , 
smh z smh z 

where coshz — A. Finally, the transmission can be re- 
lated to the evolution matrix as 



TNiE,ky) 



li^% 



(5) 



We are now in the position to discern whether an elec- 
tronic mode is characterized by an oscillating or evanes- 
cent behavior far from the barrier superlattice. To this 
aim, we preliminarly observe that, depending on E and 




FIG. 2: (Color online) Single electron transmission Ti(E, ky), 
Eq. ((5]) across a single barrier (A'^ = 1, ^ = 25 nm), as a 
function of scaled transverse wave vector hvpky/Ei and scaled 
energy E/Ei, Eq. ((BJ, with Ei « 40 meV. Here, strain is 
applied along the armchair direction, 9 = 7r/2, and we set 
£_ = 0.02, e+ = 0, and U± = 0. Cyan dashed lines delimit 
cones corresponding to the (deformed) Dirac cones outside 
(left cone) and within (right cone) the barrier (regions I-I-III 
and II, respectively, in Fig. [l]). Solid lines outside the left 
Dirac cone correspond to bound modes. 




FIG. 3: (Color online) Single electron transmission T^>(E, ky), 
Eq. ((5]) across a superlattice of A'' = 5 identical barri- 
ers (Fig. [l|, as a function of scaled transverse wave vector 
hvFky/Ei and scaled energy E/Ei, Eq. ((6]). All other pa- 
rameters are as in Fig. [2] Red lines outside the right cone 
correspond to resonant modes. 



ky, one has a propagating (resp., evanescent) wave for 
gj. < {resp., q\ > 0), where the subscript -I- refers to 
the region between two consecutive barriers [(2m — 1)^ < 
X < 2m£, m = 1,...N], and the subscript — refers to 
the region within a barrier [2(m — 1)£ < x < {2m — 1)£, 
m^l,...N] (Fig. [J). 

Fig.[2]shows the single electron transmission Tn{E, ky) 
across a single barrier, Eq. ([5]) with A^ = 1, as a func- 
tion of the transverse wave vector hvpky/ Ei and energy 
E/Ei, each scaled by the characteristic energy 
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nhvp 
2^7 '' 



(6) 



where 7= |[(1— Aa;£+)~"^+(l— Aa;£:_)~^]. Here and in the 
foUowing, strain is apphed along the armchair direction, 
9 = 7r/2, and we set e_ = 0.02, e+ = 0, and U± — 0. 
In Fig. [2J cyan dashed hnes delimit the two (deformed) 
Dirac cones defined by q^ < (left cone) and gi < 
(right cone), corresponding to regions I+III and II in 
Fig. H]), respectively. One finds that Ti{E,ky) is defined 
within the left cone and is exponentially vanishing within 
the intersection between both cones. This corresponds 
to having propagating modes in all the three regions. In 
this case, resonant modes, i.e. propagating modes with 
unit transmission, are characterized by the condition for 
stationary waves 



rmr, 



(7) 



where q- = ig_, and m is an integer. 

Fig.[3]shows the single electron transmission Tm{E^ ky) 
across a superlattice composed of five identical barriers, 
Eq. ([5]) with N = 5. Again, nonzero values of the trans- 
mission are to be found within the intersection of the 
Dirac cones corresponding to the region inside a barrier 
and between two consecutive barriers. However, at vari- 
ance of the case iV = 1, because of multiple scatterings, 
a nonzero transmission is also possible beyond the cone 
q^ < 0. This corresponds to having evanescent modes 
within the barriers. Such a phenomenon is analogous to 
what happens to photons propagating across a ID pho- 
tonic crystal with alternative layers of a left-handed and 
a right-handed material (ID metamaterial)^. As for res- 
onant modes, TN{E,ky) = 1, besides the ones given by 
Eq. d?]) regardless of N, additional resonant modes are 
given by the condition 



A = cos — 

N 



j = l,...N~l, 



(8) 



where A is defined by Eq. (j2b|) . and |A| < 1. The latter 
condition implies that these resonant modes have glob- 
ally propagating behavior. Comparing Figs. [2] and [3l one 
finds that, in the domain within both Dirac cones, in 
addition to the resonant modes given by Eq. ([7]), in the 
case A'^ > 1 there exist A^ — 1 new resonant modes given 



by Eq. ([8]). It should also be noted that in the domain 
within the left cone but outside the second the resonant 
modes, which are only given by Eq. ([8]), are characterized 
by quite narrow lines in the transmission plots. 

Outside the left Dirac cone, it is still possible to find 
bound state&i^i21iSS. Within the transfer matrix method, 
these are given by the condition^S [M^^)]^^ = 0. For 
q^ > one finds evanescent modes outside the barriers, 
and therefore also far from the superlattice structure. In 
the case A^ = 1, one finds several such confined modes 
within the second cone (Fig. [21 solid lines outside the left 
cone), whereof only one such mode survives in the region 
outside both cones. The latter is the surface mode ana- 
lyzed in Ref. [l^ In the case A^ > 1 (Fig. |3l solid lines 
outside the right cone), one finds that to each bound 
mode in the single barrier case there correspond exactly 
A^ bound states outside either cones. This is remindful 
of electron bands in solids, where the overlap of A^ peri- 
odically arranged atomic orbitals give rise to a band of 
N states. 

In conclusion, we have found that a strain-induced su- 
perlattice in graphene can accomodate additional reso- 
nant quasiparticle states, analytically characterized by 
Eq. (|5]), besides the ones usually found across a single 
barrier, given by Eq. ([7]). One finds that applied strain 
modifies the kinetic part of the quasiparticle Hamilto- 
nian, which preserves its Dirac-like character, but around 
shifted and deformed Dirac cones. This can be described 
in terms of a coordinate-dependent, periodic, profile of 
the Fermi velocity, which produces coherent effects on the 
quasiparticle transmission. Specifically, we find resonant 
modes with globally propagating behavior far from the 
superlattice, for conserved energy and transverse momen- 
tum within the intersection of the two deformed Dirac 
cones corresponding to the two alternating strained re- 
gions. Other modes are exponentially suppressed, and 
we also discuss the spectrum of bound states, which ar- 
range themselves as 'bands', depending on the overall 
number of barriers making up the superlattice. We thus 
surmise that a strain-induced superlattice in graphene 
can be used as a filter for the resonant modes here dis- 
cussed. 
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